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AMEND^^ENTS TO THE SPEClFICATtONS 
Please amend paragraph (0005] as indicated: 

[00051 A conventional way to determine a localized specinim is to use the method of fia 
Short Time Fourier Transform (STFT). In STFT, tapered moving windows of the time 
domain signal are used to compute their Fourier spectra. The general foim of this 
trajififonn i$: 

-00 

where/(l5) is the time-domain function, g(t) is the window function^ and e"^^ the 
Fourier kernel. This method has heen summarized hyNawab and Qauatieri (see Nawab 
et al-, 1988. **Short-time Fourier transform" in Advanced Topics in Signal Processiiig, 
Lim» If and Oppenheim^ Al. Eds.: Prentice Hall Signal Proce$$ing Series, 289-337.) and 
employed in practice by Par^a et al. (see Partyka et al., 1999. ^Inteipretational aspects 
of spectral decomposition in reservoir characterization: The Leading Edge, IS, 3S3-360.) 
and Marfurt et al. (see Marfurt et al., 200L "Nairow-band spectral analysis and tfain-bed 
tuning: Geophysics, 66, 1274-1283.) The analysis window function plays an important 
role in the STFT> where resolution issues are dependent to a large degree on window size. 
The longer the window size is in time^ the better the resolution of the local spectrum in 
the frequ^cy domain, but the worse the resolution in the time domain. The shorter the 
window dze is, the better the resolution of the local spectrum in the time domain, but the 
worse the resolution in the frequency domain. 
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Please amend paragraph [0007] as indicated: 
(0007] In one method, for an observed sequence 



M 

where aM,m is the AR parameter m of the A<^* order AR process, represents the nos 
components. In terms of vectors, 

where 

and rdenotes the vector transpose. If each side of the above equation is pr&-multiplied 
by complex conjugate vector X^^^ and the expected value taken, then 



where 



= (Af + l)(Af+l) 
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is the Tocplitz autoconclation matrix, ^ ■'' ' is the autocorrelation functioaat 



lag time 
density. 



ime t. f '^^ and '^[^'"^" ] Js the white noise power spottal 
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Please amend paragraph [0008] as indicated: 

[0008] By definition, the white noise is uncorrelatcd with all for m<n. Thus 

It is proven that the power spectral density Sx(f) of the series is related to the input noise 
pow^ spectral density by 



Pm 



M 

1 + E^Ji/^ exp(-z2^/mA/) 
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Please amend paragraph (0010] as Indicated: 

[0010] The latest approach to precluding or diminishing the problems due to windowing 
involves the use of the wavelet analysis* Wavelet analysis is a newly established (since 
the late 1980s) field in mathematics and signal processing. Like the Fourier transfonn, 
the wavelet transform also convolves through a discrete summation or continuous 
integration the time function (signal) with a kernel function. A method of using wavelets 
is found in Chakraborty et al, (sec Chakraborty et al., 1995, Frequency-time 
decomposition of seismic data using wavelet-based methods: Geophysics, 60, 1906- 
1 916.) as well as Xia (Xia, L., 1999, Spectral analysis of seismic data using wavelet 
transform: M.S. Thesis, University of Oklahoma.). In wavelet analysis, a wavelet is used 
as the kernel function in place of the Fourier kernel. For example, given a fimction/(1r^> 
its Fourier transform is: 

F(f)=]f(t)e"^dt (2) 

wh^e and is the Fourier kernel. Its carresponding wavelet transfomi 

is: 



where is the complex conjugate of wavelet Wavelet ^ is a function that is squaie- 
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integrable ^ ^ having zero mean and localized in both time and frequcncj, 



Therefore, represents a family of wavelets that satisfies the conditions: 



where and . aisiefeEiedtoasthescaleofthewaveletandrisiQfensdto 



as the transLation param^er. Note that the wavel^ is nonnalized such that ||^|| = 1 . 
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Please amend paragraph [0012] as indicated: 



[0012] In MPD, a family of wavelets is defined by the forai 



where 



and 



. a is referred to as the scale, r a$ the translation parameter, 



and ^ as the frequency modulation. Each wavelet in the family is called a time-jfrequeacy 
atom. If ^(t) is Gaussian, these atoms are called Gabor atoms. As shown by Mallat et aL, 
Gabor atoms provide excellent time-frequency resolution. These atoms have 
combinations of all possible time and frequency widths and as a result constitute a 
redundant set. Once atoms are defined^ a best match between the signal and these atoms 
is found by projecting the atoms onto the signal and then computing the maximum. A 
residue is then computed by subtracting from the original signal the product of the atoms 
and the cross product of the selected atom and the signal. This decomposition is 
continued until the energy of the residue falls below some ttuneshold. This method has the 
ability to obtain a good resolution in both time and frequency for data within an 
mtcncncdiatc frequency range like seismic data. 
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